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We study the routing of quantum information in parallel on multi-dimensional networks of tunable 
qubits and oscillators. These theoretical models are inspired by recent experiments in superconduct- 
ing circuits using Josephson junctions and resonators. We show that perfect parallel state transfer is 
possible for certain networks of harmonic oscillator modes. We further extend this to the distribution 
of entanglement between every pair of nodes in the network, finding that the routing efficiency of 
hypercube networks is both optimal and robust in the presence of dissipation and finite bandwidth. 
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The rapid development of coherent superconducting 
quantum bits (qubits) [1 demands continued theoretical 
analysis of quantum networks. Current experiments [2 
are in fact realizing multi-dimensional coupling topolo- 
gies, and advanced coupling networks are on the horizon 
HSU- While these can provide the hardware to run quan- 
tum algorithms, they can also be seen as the fabric of a 
quantum network to route quantum information between 
arbitrary nodes in parallel. This network can be used for 
the transfer of entanglement, which can be locally pu- 
rified [5 for teleportation, error correction [6 or other 
tasks in quantum information processing. 

Efficient routing of entanglement between nodes can be 
be accomplished provided two conditions are met. First, 
it must be possible to program the network to transfer 
quantum states between arbitrary nodes. This is possi- 
ble if the qubit frequencies or inter-qubit couplings can 
be dynamically controlled. Second, the network should 
exhibit high fidelity state transfer between nodes in par- 
allel, as a multi-user network. This is in contrast to the 
various methods studied for state transfer in a single-user 
network; see [7] for a review. Here we consider two such 
networks that meet these conditions — the hypercube and 
the complete graph, illustrated in Fig. 1. 

To exploit the full parallelism of these networks, we 
introduce a new approach, using oscillator networks to 
route multiple excitations between nodes. Previous the- 
oretical work on quantum routing [8 used single state 
transfer of qubits, while the study of entanglement dy- 
namics in oscillator networks [9] has commonly focused 
on single-excitation or continuous- variable entanglement 
(see [10] for a notable exception). 

In this Letter we theoretically analyze parallel state 
transfer and quantum routing in these networks. We be- 
gin by providing the first calculations for the fidelity of 
parallel state transfer, showing that multiple nodes can 
faithfully send entanglement through the network at the 
same time. Using these results, we consider entangle- 
ment distribution schemes to evaluate the routing effi- 
ciency, and provide analytical results for the complete 
and hypercube networks in the presence of cross-talk and 
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FIG. 1: Parallel state transfer on programmable quantum net- 
works. Each node is an oscillator with a tunable frequency. 
Each line (solid or dashed) indicates a coupling between os- 
cillators. Solid lines indicate couplings between oscillators 
with the same frequency; dashed lines indicate couplings be- 
tween oscillators with different frequencies. High fidelity state 
transfer occurs for large detuning, (a) Hypercube network 
with d = 3, programmed into two subcubes (red and blue 
squares). Each node is labeled by a bit-string of length d = 3, 
here with the first m — 1 bits indicating the subcube. In the 
qubit-compatible scheme (QC), one entangled pair is sent on 
each subcube, as indicated by the arrow for the inner (red) 
square, (b) In the massively parallel scheme (MP) scheme, 
multiple entangled pairs are sent between every node of each 
subcube, as indicated by the arrows for the inner (red) square, 
(b) Completely connected network with N = 8, programmed 
into N/2 = 4 two-site networks. 



dissipation. This lays the foundation for analyzing many 
future quantum networks, and can guide future experi- 
ments with superconducting or other qubit-oscillator sys- 
tems. 

General framework. A quantum network is described 
as a graph Q = (V,E) of vertices V connected by edges 
E, with a Hamiltonian T~L given by 

H = h^2 Qv V aia v + ft ^2 ^uv {o) u a v + a\a u ) . (1) 
vev {u,v}eE 

where a\ and a v are the creation and annihilation op- 
erators for the oscillator at vertex v. We confine our 
attention to the case when the couplings between nodes 
are all equal to Qq, although we do require that frequen- 
cies of each node Q vv be programmable from one use of 
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the network to the next. Such networks can be realized 
by tunable superconducting qubits coupled by capacitors 
or through resonators, as realized in recent experiments 
[2]. While tunable couplings have recently been demon- 
strated [11] , their use for quantum routing requires addi- 
tional study. Quantum oscillator networks have not yet 
been realized, but great progress has been made towards 
qubit-oscillator systems [12] and tunable oscillators [T3] . 

To analyze the transfer of entanglement between nodes 
of a quantum network, we consider the following modi- 
fication of the state transfer protocol [7J [14]. At each 
node of this network there is a register of auxiliary qubits 
capable of generating pairs of qubits in the Bell state 
|$+) = (|00) + \ll))/y/2. Some subset of the nodes 
{si . . . %} act as senders, and some subset {r± . . . tm} 
act as receivers. Each sender prepares his auxiliary qubits 
in a Bell state, and then swaps one of these qubits onto 
the network (as done in recent experiments [12]). The 
goal of sender Sj is to faithfully transfer half of her en- 
tangled pair to receiver rj. 

More formally, we consider the time-evolution of the 
initial state 

M 

\9{t = 0)> = 2- M / 2 [J (l + at .&J .) | vac ), (2) 
i=i 

where b\ is the raising operator for the auxiliary qubit at 
vertex v and the vacuum state |vac) is the state where 
all oscillators and qubits are in the ground state. This 
initial state is a product of Bell states for each sender's 
qubit-oscillator pair. We evolve this state for some time 
T, |tf (T)) = E/(T)|tf (0)) with U(T) = e- iUT l n : 

M 

|*(D> = 2- M / 2 n (l + U{T)a\ 3 U\T)b\^ |vac). (3) 

Our target state has the entanglement transferred to the 
receiver's node: 

M 

^target) = 2~ M '* JJ (l + 4$,) \™) > (*) 
3 = 1 

where e l ^ j is a known phase that can be corrected by local 
unitary operations on the pair {sj,rj). To characterize 
the state transfer, we use the fidelity of the final density 
matrix pj shared between receiver rj and the sender Sj 
(after the appropriate local correction) with our intended 
perfect Bell state: 

F Sj ^ rj = ($ + | ft |* + >. (5) 

We calculate a lower bound on these fidelities by ana- 
lyzing the time-evolved state and taking its overlap with 
an appropriately chosen product state |^ ) = 
This calculation is significantly simplified by the linearity 
of the time-evolution of the oscillator modes, for which 

U(t)a J U\t) = ^K jk (t)a k , (6) 

k 



where we have defined the unitary mode evolution matrix 
K(t) = exp(— iQt) and where the N x TV coupling matrix 

has elements fl uv . Note that perfect state transfer 
would occur for U(T)a\.W(T) = e^a\ , 

Parallel state transfer on the hypercube. Christandl et 
al. showed that one can perform perfect state transfer 
from corner-to-corner of a d-dimensional hypercube in 
constant time T = 7t/(2Qq) [15]. Here we analyze how 
this result can be extended to the transfer of quantum 
states in parallel, by splitting the cube into subcubes. 
Specifically, by tuning the frequencies of each node, the 
d-dimensional hypercube can be broken up into 2 m sub- 
cubes each of dimension d — m, as shown in Fig. [TJa) 
for d = 3 and m = 1. These subcubes can be made 
to act as good channels between their antipodal nodes 
by separating the oscillator frequencies for each channel 
from adjacent channels by an amount Auj. For fixed cou- 
plings, there is still the potential for cross-talk between 
channels, which we now analyze. 

Perfect state transfer on the hypercube is most conve- 
niently described using a binary labeling scheme [3], in 
which each node is labeled by a bit-string of length d. For 
parallel state transfer, this string separates into m bits 
that specify the various subcube channels, and (d — m) 
bits for the position on each subcube, as illustrated in 
FigfTJa). Using such a labeling, we choose the detunings 
so that the coupling matrix takes the form 




where and are Pauli matrices for the 2 d nodes 
of the graph. This coupling matrix has the nice property 
that all oscillators on the same subcube have the same 
frequency, and are detuned from adjacent subcubes by 
Auj. Crucially, we have written Eq. ^ as a sum of d 
commuting matrices. For Auj ^> Qq, the first m terms of 
ft keep the subcubes separate, while the remaining (d — 
m) terms serve to transfer excitations between corners 
on each subcube. 

For this choice of fJ, the mode evolution matrix K(t) = 
exp(— iVtt) can be evaluated exactly. Then, by choosing 
the M = 2 m senders and receivers to transfer in the same 
direction along each of their subcubes (e.g. from 00 • • • 
to 11 • • • 1), one can show [16] that the fidelity satisfies 

F.^r, >l-lmr] 2 sin 2 + 0( V 3 ), (8) 

for all j, where 77 = 2Qq/Auj and £t = \A + V~ 2 ' Note 
that for large detuning (77 <C 1) each sender transfers 
entanglement to her corresponding receiver with high fi- 
delity. In addition, the fidelity exhibits resonances due 
to the dependence on sin 2 by picking 77 just right, one 
can make sin 2 £p = 0. Remarkably, the exact calculation 
shows that for these detunings the fidelity is equal to one 
to all orders in 77, realizing perfect parallel state transfer. 
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FIG. 2: The fidelity of entanglement transfer on the hyper- 
cube as a function of the detuning parameter n = 2Qq/Alu. 
The qubit curves are, from top to bottom (for small 77), nu- 
merical simulations for dimension d = 2 — » 6. Each network 
is split into M — 2 subcube channels with one sender and re- 
ceiver per channel, with entanglement being sent in the same 
direction on each channel. Also shown also is the lower bound 
of Eq. (|8| for the oscillator network with m = 1. 

Equation ([8| was derived for entanglement transfer on 
oscillator networks. However, as long as only one sender- 
receiver pair uses each channel at a time, numerical cal- 
culations, shown in Fig. [2j indicate that qubit networks 
behave similarly. For this reason we call the parallel state 
transfer protocol discussed so far the "qubit-compatible'' 
(QC) protocol. There are some notable differences be- 
tween qubits and oscillators, namely qubit s do better on 
average, but do not exhibit perfect state transfer. 

Besides the ability to perform truly perfect parallel 
state transfer, oscillator networks have another feature 
that qubit networks lack: the capacity for massively- 
parallel (MP) entanglement transfer. Because multiple 
excitations on an oscillator network pass through each 
other without interacting, every node on an oscillator 
hypercube network can act as both sender and receiver 
at once, with each transferred pair having a fidelity given 
by Eq. (J8|. That is, instead of sending M = 2 m states in 
parallel (as described above), one can send 2 d states all 
at once, by transferring 2 d ~ m states on each of the 2 m 
subcubes. This is illustrated by the arrows in Fig. [ljb) 
for d = 3 and m = 1, with 4 states transferred on each 
square. In this scheme, all of the network is used all of 
the time. The fidelity of each transfer is precisely the 
same as in the QC scheme. As we will soon show, this 
allows for optimal routing efficiency on the hypercube. 

Parallel state transfer on the complete graph. The com- 
plete graph is intrinsically qubit-compatible and mas- 
sively parallel, as all pairs are simultaneously coupled. 
Such a network can be realized by coupling many su- 
perconducting qubits to a common metallic island or a 
high-frequency superconducting resonator [2]. One can 



address different pairs of qubits by tuning them into res- 
onance, and separated from other pairs by Auo. Using 
perturbation theory on such a configuration, one finds 
that 

7T 2 

^complete >l-y^ + 0(^ 3 ), (9) 

where again rj = 2Qq/Auj. In contrast with the hy- 
percube, here there is no dependence on the number of 
nodes, and the optimal transfer time is in fact different 
for different pairs. The latter complication is ignored in 
the following. However, the dependence on the band- 
width through 77 plays an important role in the routing 
efficiency. 

Entanglement Distribution and Routing Efficiency. 
Using parallel state transfer, these networks can be pro- 
grammed for the efficient transfer of quantum informa- 
tion. To quantify this efficiency, we consider the distri- 
bution of entanglement between every node. That is, we 
consider protocols that route entanglement between ev- 
ery pair of nodes on the network. In addition, we consider 
the realistic condition that all oscillator frequencies are 
within a finite bandwidth. We now calculate the routing 
efficiency for our three protocols. 

To quantify the routing efficiency, we use the rate of 
entanglement distribution. We define the distribution 
time Tr> as the time it takes for all nodes to share an 
approximate Bell pair, and the distribution rate as the 
total number of entangled pairs shared, weighted by their 
fidelities, divided by the distribution time: 

n = ^ E ( 10 ) 

pairs {u, v} 

This rate distinguishes one-dimensional schemes [14], 
that take a long time to transfer a small amount of en- 
tanglement, from multi-dimensional schemes that take a 
short time with perfect state transfer [15]. Note that all 
serial schemes with transfer time T have 7Z < 1/T. By 
using parallel state transfer, as described above, this rate 
can be made even greater. 

For the complete graph, we split the network into N/2 
pairs. Each pair is detuned by Auj from adjacent pairs, 
for a total bandwidth of co> max — c<j m i n = NAuj/2. The 
distribution time Tp is simply the swap time T = 7r/2^ 
times the number of steps required to couple each pair of 
nodes, or Tjj — T(N — 1). The total amount of entan- 
glement shared is N(N — 1), or two for each pair. Using 
the fidelity of Eq. ([9| , we find that 

For the hypercube, we split the network into 2 m sub- 
cubes. Here only adjacent subcubes need be detuned, 
for a total bandwidth of cj max — co> m in = mAcj, as in- 
dicated in Eq. ([7]). To calculate the distribution time, 
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we consider our two schemes separately. For the qubit- 
compatible scheme, in which one excitation is on each 
subcube channel, the split network must be used for ev- 
ery set of antipodal pairs, or 2 d ~ m ~ 1 times. Given that 
there are (^) such splittings, and that each transfer takes 
time T, the total distribution time is 



m=0 v 7 



)d—rn — l 



= -(3 d -l)T. (12) 



Using Eq. M with sin £ T = 1, the sum in Eq. (10) can 



be performed analytically, from which we find, for large 
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d 2 (d + 3) , 



(13) 

where N = 2 d and the prefactor is (4/3) d = N lo ^ 3 \ 

For the massively parallel scheme, every antipodal pair 
transfers simultaneously in time T, giving 



P (MP) 
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with a rate 
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(2 d - 1)T, 



(14) 
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;d 2 (d + 3) . (15) 



These rates are conservative estimates; higher rates are 
possible by using the resonances seen in Fig. [2] 

These three distribution rates are plotted as a function 
of the number of nodes in Fig. [3j where we have fixed the 
bandwidth appropriate to recent superconducting qubit 
experiments [2]. For the hypercube schemes, the mas- 
sively parallel protocol is more than quadratically better 
than the qubit-compatible scheme. Entanglement trans- 
fer on the complete graph quickly fails due to significant 
cross-talk for N « 20. One might expect that this is due 
to the large number of couplings, but it is actually due 
to the finite bandwidth of the network. It is clear that 
studying extended coupling schemes such as the cavity 
grid [4] is an important task. 

Decoherence and Disorder for the Hypercube. Experi- 
mental issues related to hypercube state transfer, includ- 
ing decoherence and disorder, have been analyzed previ- 
ously [3]. These results can be applied directly to the QC 
scheme. Decoherence will simply reduce the fidelity (and 
1Z) by a factor ~ e~ T l Tci for arbitrarily large subcubes 
[3], where T2 is the total dephasing time. A full analy- 
sis of decoherence for the multiple excitations of the MP 
scheme will require additional study, but the dominant 
source of decoherence for superconducting oscillators is 
dissipation [17]. For this decoherence process, one can 
calculate that the fidelity will be reduced by a factor 
~ e _T / Tl , where T\ is the dissipation time. 



10 4 



1 1000 

a 

I 100 



10 



MP Hypercube 



QC Hypercube 



Complete Graph 



10 100 1000 10 4 

Number of Nodes 



FIG. 3: Entanglement distribution rate 1Z (in units of 1/T) 
as function of the number of nodes N in a quantum network. 
Three distribution schemes are shown: the massively paral- 
lel (MP) and qubit-compatible schemes on the hypercube of 
dimension d (each with N — 2 d ), and the complete graph 
of size N. Each network was chosen to to have a coupling 
of Qo/2n — 20MHz with a bandwidth (cj max — cj m i n )/27r — 
2GHz. 



In conclusion, we have analyzed how entanglement can 
be routed between all nodes of the hypercube and com- 
pletely connected networks. This has been accomplished 
by parallel state transfer and analytical calculations of 
the entanglement distribution rate. In the ideal case, 
oscillators on both the hypercube and complete graphs 
achieve optimal efficiency. This efficiency is robust for 
the hypercube in the presence of finite bandwidth and 
dissipation. These results provide further evidence that 
superconducting resonators are an important element for 
quantum information processing, and motivate further 
study of parallelism in quantum networks. 
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